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Abstract. The large-n expansion is developed for the study of critical behaviour 
of d-dimensional systems at m-axial Lifshitz points with an arbitrary number m 
of modulation axes. The leading non-trivial contributions of 0(l/n) are derived 
for the two independent correlation exponents n^2 and «i4i and the related 
anisotropy index 8. The series coefficients of these 1/n corrections are given 
for general values of m and d with < m < d and 2 + m/2 < d < 4 + m/2 in 
the form of integrals. For special values of m and d such as (m, d) = (1, 4), they 
can be computed analytically, but in general their evaluation requires numerical 
means. The 1/n corrections are shown to reduce in the appropriate limits 
to those of known large-n expansions for the case of d-dimensional isotropic 
Lifshitz points and critical points, respectively, and to be in conformity with 
available dimensionality expansions about the upper and lower critical dimensions. 
Numerical results for the 1/n coefficients of r)^, VLi an d 8 are presented for the 
physically interesting case of a uniaxial Lifshitz point in three dimensions, as well 
as for some other choices of m and d. A universal coefficient associated with the 
energy-density pair correlation function is calculated to leading order in 1/n for 
general values of m and d. 
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1. Introduction 

Systems exhibiting critical behaviour can be divided into universality classes such that 
all members of a given class have the same universal critical properties (Fisher 1983). 
The universality classes are represented by field theories, such as the n-component </> 4 
models in d space dimensions, which are nontrivial whenever d is less than the upper 
critical dimension d* above which Landau theory holds. Proper analyses of such field 
theories usually require sophisticated tools such as renormalization group approaches, 
elaborate numerical simulations or a combination of both (Wilson and Kogut 1974, 
Domb and Green 1976). 

The commonly employed successful analytical methods are dimensionality 
expansions about the upper and lower critical dimensions d* and di (Wilson and 
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Fisher 1972, Polyakov 1975, Brczin and Zinn-Justin 1976a, 1976b, Bardeen et al 
1976), the massive renormalization group (RG) approach in fixed dimensions (Parisi 
1980), and the 1/n expansion in inverse powers of the number n of order-parameter 
components (Ma 1973, Abe 1973). Remarkably precise estimates of critical indices 
and other universal quantities of three-dimensional systems have been obtained both 
via e = d* — d expansions and the massive RG approach at fixed d by computing 
perturbation series to sufficiently high orders and resumming them (Guida and Zinn- 
Justin 1998, Pelissetto and Vicari 2002, Privman et al 1991). Expansions about the 
lower critical dimension dg , on the other hand, seem to have a more modest potential 
for precise estimates, unless they are combined with information from other sources. 

In this paper we will be concerned with the 1/n expansion. Our aim is to develop 
this approach for the study of critical behaviour at m-axial Lifshitz points. A Lifshitz 
point (LP) is a multicritical point at which a disordered phase meets both a spatially 
uniform ordered and a spatially modulated ordered phase (Hornreich et al 1975a, 
Hornreich 1980, Selke 1992, Diehl 2002, 2004). The disordered phase is separated from 
the two ordered ones by a line of critical temperatures T c = T c (g), depending on a 
thermodynamic nonordering field <?, such as pressure or a ratio of an antiferromagnetic 
and a ferromagnetic coupling. The LP, located at Tl = T c (gL), divides this critical line 
into two sections. In the modulated ordered phase, the wave-vector associated with 
the modulation, q mo d = Qmod(T, g), varies with g and temperature T. Without loss of 
generality, we can consider the ferromagnetic case. Then uniform order corresponds 
to q mo d = 0, and q mo d(T,g) vanishes at the LP. The LP is called m-axial if the 
wave- vector instability that sets in at the LP occurs in an m-dimensional subspace of 
d-dimcnsional space, i.e. q mo d € K m with < m < d. The limiting values m = and 
m = d correspond to the cases of a usual critical point (CP) and the isotropic LP, 
respectively. 

As is well known, the large-n expansion leads to self-consistent equations (Abe 
1973, Brezin et al 1976, Vasiliev et al 1981a). These are considerably more difficult to 
handle than ordinary perturbation expansions, which is the essential input required 
in the alternative approaches mentioned above (dimensionality expansions, massive 
RG approach). For this reason the available series expansions in 1/n are restricted 
to low orders, even for the much simpler case of a conventional CP. Furthermore, 
1/n expansions normally converge slowly. || On the other hand, this technique has 
a number of very attractive features. One is its capability of treating fluctuation 
effects in a systematic, nonperturbative manner. % Another is that it can be applied 
for any fixed value of d between dt and d*. No additional expansion in a small 
parameter such as e = d* — d is required. Owing to these appealing properties the 
1/n expansion continues to be an important tool for the analysis of nontrivial field 
theories. The spectrum of problems to which it has recently been applied and made 
significant contributions is impressively rich. It ranges from the study of classical 
spin models (Brczin 1993, Vasiliev 1998, Zinn-Justin 1996, Moshe and Zinn-Justin 
2003, Campostrini et al 1998, Pelissetto et al 2001, Gracey 2002a, 2002b), critical 
behaviour in bounded systems (McAvity and Osborn 1995), the physics of disordered 

|| According to some recent results (Baym et al 2000, Arnold and Tomasik 2000), the shift of the 
critical temperature of a dilute Bose gas seems to be an exception to this rule. 

1 A familiar numerical RG scheme having this capability employs the so-called "exact" RG equations 
(Wilson and Kogut 1974). Its only application to critical behaviour at LPs we are aware of is the 
recent work by Bcrvillicr (2004). This deals with the case of a uniaxial LP with n = 1 in d = 3 
dimensions and uses the lowest ("local potential") approximation. 



Large-n expansion. 



3 



elastic media (Le Doussal and Wiese 2003, 2004) and models of stochastic surface 
growth (Doherty et al 1994) to problems of high-energy physics (Aharony et al 2000, 
Moshe and Zinn- Justin 2003) and quantum phase transitions (Franz et al 2002). 

Unfortunately for critical behaviour at LPs, hardly any results have yet been 
worked out by means of it. Previous large-n analyses for general values of m G 
[0, d] (Kalok and Obermair 1976, Hornreich et al 1977, Mukamel and Luban 1978, 
Frachebourg and Henkel 1993) were restricted to the limit n — > oo. To our knowledge, 
the only exceptions in which corrections to order 1 /n were computed deal exclusively 
with the case m = d of a d-axial LP (Hornreich et al 1975b, Nicoll et al 1976, Inayat- 
Hussain and Buckingham 1990). 

This lack of results for general m is due to the complications arising from the 
combination of two special features such systems have: (i) the anisotropic nature of 
scale invariance they exhibit at LPs, and (ii) the unusually complicated form of the 
two-point scaling functions they have in position space already at the level of Landau 
theory (Diehl and Shpot 2000, Shpot and Diehl 2001). Anisotropic scale invariance 
means that coordinate separations dz a along some directions scale as nontrivial powers 
of their counterparts drp along the remaining orthogonal ones, i.e. dz a ~ [drpf , where 
6, called anisotropy exponent, differs from 1. Of course, anisotropic scale invariance 
is a feature encountered also in studies of dynamical critical behaviour at usual CPs, 
where time t scales as a nontrivial power of the spatial separations, and for which 
analyses to 0(l/n) can be found in the literature (see e.g. Halpcrin et al (1972); for 
more recent uses of the 1/n expansion in dynamics, see e.g. Bray (2002) or Moshe and 
Zinn- Justin (2003)). It is the combination with (ii) that makes consistent treatments 
of fluctuation effects so hard for critical behaviour at to- axial LPs. The very same 
difficulties had prevented the performance of full two-loop RG analyses within the 
framework of the e = d* — d expansion for decades, and produced controversial 0(e 2 ) 
results (Mukamel 1977, Sak and Grest 1978, Hornreich and Bruce 1978, Hornreich 
1980, Mergulhao and Carneiro 1999) until such a RG analysis for general m was 
finally accomplished in (Diehl and Shpot 2000, Shpot and Diehl 2001).+ 

In order to overcome the mentioned technical difficulties, we shall adapt and 
extend the elegant technique suggested by Vasiliev et al (1981a) and utilized in 
subsequent work (Vasiliev et al 1981b) to compute the 1/n expansion of the standard 
critical exponents r] and v for the case m = of a CP to 0(l/n 2 ). (In Vasiliev et al 
(1981c), the authors managed to compute rj even to 0(l/n 3 ) by means of a conformal 
bootstrap method.) We proceed as follows: 

We start in section 2 by specifying the expected scaling forms of the required two- 
point functions. These are employed in section 3 to solve the resulting self-consistent 
equations to the appropriate order 1/n. Matching the anticipated asymptotic large- 
distance behaviour of the cumulants with the solutions of these equations yields 
consistency conditions. From these, the 0(l/n) terms of the two independent 
correlation exponents r\Li and t\la can be determined for general values of to in the 
form of finite, numerically computable integrals. 

In section 4 we verify that these coefficients reduce in the limits to — > and 
m — > d to the analytical expressions for the CP case and the isotropic LP obtained 
in (Ma 1973, Abe and Hikami 1973) and (Hornreich et al 1975b), respectively. In 
section 5 we consider the behaviour of our 0(l/n) results as e = A + m/2 — d^> 0, and 

+ For a discussion and critical assessment of recent work (dc Albuquerque and Lcitc 2001, Leite 2003) 
giving alternative results, see (Diehl and Shpot 2001, 2003). 
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explicitly reproduce the large-n limits of the 0(e 2 ) coefficients of (Diehl and Shpot 
2000, Shpot and Diehl 2001) for the correlation exponents 77^2 and 77^4 in the case 
m = 2 of a biaxial LP. The consistency of our results with known expansions about 
the lower critical dimension de — 2 + to/2 is demonstrated in section 6. 

In section 7 we first consider the special cases (m,d) — (1,4) and (4,5). For the 
former, we derive fully analytical expressions for the 1/n coefficients of the correlation 
exponents and the anisotropy exponent. The latter case requires some numerical work, 
even though part of the calculations can also be done analytically. Next, we turn to 
the physically interesting case (to, d) — (1, 3) of a uniaxial LP in three dimensions and 
present numerical results for the coefficients of the 1/n contributions. In section 8 we 
compute the universal scaling function of the energy-density pair correlation function 
and a related universal amplitude to leading order in 1/n. We give analytical results 
for the expansion coefficient of the latter for general values of (to, d) with de < d < d* . 
Our results are briefly discussed and put in perspective in the closing section 9. Finally, 
there are two appendixes describing technical details. 

2. Scaling properties of the two-point functions 

We wish to consider the theory of an n-component vector field <p{x) = 
(<pi(x), . . . , <p n (x)) in d space dimensions. Writing x = (r,z), we decompose the 
position vector x G R d into a "parallel" component z e M m and a "perpendicular" 
one r € R m , where to = d — to is the codimension of to. Likewise, d r = d/dr 
and d z = d/dz denote the corresponding components of the gradient operator 
V = (d r , d z ). The Hamiltonian of the model we will investigate is given by 

H[<f>] =Wo[0]+W int [0] , (1) 

with the 0(n)-symmetric free and interacting parts 

n [<P} = l Jd d x[{d r <pf + {dl<pf + p{d z 4>f + r^] (2) 

and 

H mt [<p] = ^j& d xtf , (3) 

respectively. 

In order for the model to have a LP, the dimension d must exceed de(m, n), the 
lower critical dimension, which is believed to be d^ n \m) = 2 + to/2 for the case 
n > 1 of continuous 0(n) symmetry we are concerned with here (Hornreich et al 
1975a, Grest and Sak 1978, Diehl 2002). 

A cautionary remark should be added here. The arguments of Hornreich et al 
(1975a) and the expansion about the dimension d = 2 + to/2 (Grest and Sak 1978) 
actually show only that the homogeneous ordered phase becomes thermodynamically 
unstable to low-energy excitations at temperatures T > whenever d < 2 + to/2. 
To establish that a LP exists for d > 2 + to/2, one would have to prove additionally 
the existence of a modulated ordered phase that is separated via a second-order line 
from the disordered phase. According to an argument by Garel and Pfeuty (1976) 
one expects the transition from the disordered to the modulated ordered phase to 
be described by a 2n-component cf> A model whose Hamiltonian has 0(n) x 0(n) 
symmetry. Although our results are consistent with the existence of a LP for n > 1 and 
d > 2 + to/2, we cannot rule out the possibility that for some values of n, fluctuations 
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might change this transition into a discontinuous one. Such a scenario apparently 
occurs in the (n = 1) case of ternary mixtures of A and B homopolymers and AB 
diblock copolymers, where fluctuations were found to transform the continuous mean- 
field transition between the disordered and lamellar phases into a discontinuous one 
(see e.g. Diichs et al (2003) and its references), in accordance with general arguments 
given by Brazovskii (1975). 

Whenever a LP exists, we denote the values of f and p at which it is located by 
tl and f>L, respectively. 

To derive the large-n expansion, it is convenient to introduce an auxiliary scalar 
field ip = ip(x) and rewrite the interaction part H.- m t[(p] as (Brczin et al 1976, Vasilicv 
et al 1981a, Vasiliev 1998, Moshc and Zinn- Justin 2003) 

e-Wmt[«] = const J VTPe-^ 6 ** [* 2 -*Va* 2 *] . ( 4 ) 
The full Hamiltonian defined by equations (l)-(3) then becomes 

H[<f>, V] =H [4>] + \f d d x [i> 2 - iVXcf, 2 V-] , (5) 

up to an irrelevant additive constant. 

In view of this reformulation of the model it is natural to consider correlation 
functions involving besides the fields <f> a with a — 1, . . . , n also tp. We need in particular 
the two-point functions. Let us consider the disordered phase. Since in it the 0{n) 
symmetry of the Hamiltonian (5) is not spontaneously broken, the mixed correlation 
functions (4> a ip) vanish. Furthermore, invariance under translations and rotations in 
the position subspaces R m and M. m implies 

(<j) ai (x + x') 4> a2 (x')) = S aia2 G (r, z) (6) 

and 

(iP(x + x')iP(x'))=G^(r,z), (7) 

where r = |r| and z = \z\ are the lengths of the perpendicular and parallel components 
of x = (r,z). Writing the wave- vector conjugate to x as k = (p,q), we introduce 
Fourier transforms G^^p^q) of the two-point functions G<p^(r,z) via 

f(d) 

G^{r,z) = / G^(p,g)e ,fc -", (8) 
Jk 

where 

f d d k 
J k = J Rd (2^p = J p J q 



(9) 



is a convenient short-hand. 

At a LP, the correlation functions are expected to decay as powers of r and 
z in the long-distance limits r, z — > oo provided df^ < d < d* . Furthermore, they 
should be anisotropically scale invariant. Thus, for appropriate choices of the values 
of the scaling dimensions and A^, of the fields 4> a and ip and of the anisotropy 
exponent 9, the limits 

lim b 2A * G (6r, b e z) = G< as) (r, z) (10) 

and 

lim b 2A * G^{br,b e z) = G { f(r lZ ) (11) 
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should exist, yielding nontrivial asymptotic functions G^(r, z). Similar results must 
hold for their Fourier transforms, namely 

Urn G$ (b~ V, b-%) = G<g (p, g), (12) 

where the scaling dimensions are given by 

2A 9 j iV , =d-m/2-2A^. (13) 
The asymptotic functions (12) are generalized homogeneous functions satisfying 
6j>, 9) = P" 2A ^ Gg(l,p-^) = q -***P (q-Wp, 1). (14) 

The exponents governing the momentum dependences of G^ are conventionally 
written as 

2A = 2-7? L2 , 2A^/9 = 4- VL4l (15) 

which defines the correlation exponents r]L2 and 77^4 • 

In the limit n — > 00, the momentum-space propagators at the LP become 

Gj )) (p,«) = (p 2 + g 4 )- 1 ) = 1, (16) 

so that the exponents take their known spherical-model values 

•lS-l2-0, *» = l/2. (17) 

For finite n, fluctuations modify both the propagators and the critical exponents. 
Making the ansatzes 

VL2 (m,d) = ^ + 0(n- 2 ), 6{m, d) = ± + — + 0(n- 2 ), (18) 
n 2 n 

we shall determine the coefficients ry£l an d 0^ of the 1/n corrections in what follows. 
For rjL4 the 1/n expansion is analogous to that of 77^2; the scaling relation (15) implies 
that its 1/n coefficient is given by 

2(^+4^)). (19) 



3. Self-consistent equations 



The self-consistent equations that must be solved in order to determine the 0(l/n) 
corrections to the correlation exponents can be derived by well-known steepest descent 
or diagrammatic methods, which are amply described in the literature (see e.g. Ma 
1973, Abe and Hikami 1973, Brezin et al 1976, Vasiliev et al 1981a, Vasiliev 1998, 
Moshc and Zinn- Justin 2003). Solving them is, however, another story: Owing to the 
anisotropic character of scale invariance at LPs, the correlation functions and G^ 
involve scaling functions, rather than being simple powers. Finding solutions to these 
equations is therefore a nontrivial task. 

We use an appropriately modified variant of the technique advocated by Vasiliev 
et al (1981a, 1981b) and employed to compute the critical exponents 77 and v of the 
standard n-component \<p\ model up to order n~ 2 . The crux of this method is to 
make a scaling ansatz for the asymptotic large-length-scale form of the solutions and 
determine the exponents it involves from the conditions implied by the self-consistent 
equations. 
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These equations can be written as 

G~J = G^-^ = G^-S >-..., (20) 



r -l _ (o(0)-l v _ ^(0)- 



■O--- (21) 



where the full and broken lines denote the full propagators (6) and (7), respectively. 

Since the propagator (6) is proportional to S aia2 , each closed loop of G^ lines 
produces a factor n. We assume that — 0(1) as n — > oo. Then it follows 
from equation (21) that G^ 1 ~ n as n — > oo. For the self-energy function E,/, an 
expansion in inverse powers of n follows, with the leading contribution being of order 
n. Consequently G^ — 0(l/n), and since G^ = 0(1), only a finite number of graphs 
contribute to E^ and E^, in calculations up to a given order in 1/n. If we restrict 
ourselves to the first nontrivial order 1/n, we may drop all contributions represented 
by the ellipsis in these two equations. 

According to equations (16), the free parts G^ _1 and G^ -1 are analytic in the 
momenta p and q. Since the asymptotic forms of the self-consistent solutions involve 
nontrivial powers of momenta for generic d, as we shall see, these free parts do not 
contribute to G^(p, q) and may be discarded. From equation (21) we thus find 

= • (22) 

while equation (20) yields 
with 

Am) Am) 

F( P ,q)= / G ( ; s \\p'+p\M' + q \)G { ; s \p',q'). (24) 

v' * Q f 

Just as G^p(p, q), F(p,q) is a generalized homogeneous function; equation (24) can 
be combined with the scaling form (14) of G^ as ' to conclude that 

F(p, q) = p ™+e™-4A, F{l lP - e q) (25) 

= q (^+0m-4A^)/e F ( q -i/e pj i) . (26) 

We can now insert these expressions together with their analogues (14) for 
G^ s ^(p,q) into equation (23), choosing cither one of the components p and q of the 
external momentum to be zero. Matching the amplitudes of the corresponding powers 
p 2A « and q 2A <t-/ e on both sides of the resulting two equations yields the consistency 
conditions 



„ Am) p 4A,-m f(m)G^(iy\p+l\-U) _ Tf . 

2 = J P \pTT^J g W^T) = /l(n) (27) 



and 



„ Am) q 4A,/0-m Am) ( ? V* \q + l|-V»p, l) 

2=1 JqTT^L fP ^ 
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where 1 designates an arbitrary unit vector in to or to dimensions. 

Since the left-hand sides are of order n, so must be the right-hand sides. Hence 
the integrals I\{n) and 7 2 (n) must have simple poles at 1/n = 0. Consider first 
I\(n). Its inner integral J^" 1 ' approaches a p- independent value as p —> oo. Recalling 
equations (15) and (18), we see that the ratio of p-dependent factors in front of 
this integral behaves as ~ p2-m+o(i/n) m ^iis limit. Had we regularized the p- 
integral by a large-p cutoff A, it would have ultraviolet (uv) divergences of the form 
i Y2-s+o(i/n)^ s _ 0,1,2. In the dimensionally regularized theory, uv divergences 
_ A 0(i/n) = i + o(l/n) lnA + ... become poles at 1/n = 0. Conversely, contributions 
from the inner integral / falling off faster than p~ 2 as n = oo do not contribute to 
the residue of the pole ~ (1/n) -1 , and hence may be dropped when calculating it. 

The upshot is that the O(n) contribution of I\{n) can be determined as follows. 
We add and subtract from G^ s \l,p 9 \p + 1\~ 9 q) its Taylor expansion in the second 
variable u p (q) = p e \p + l\~ e q about its limiting value for p — > oo, u 00 (q) = q, to 
the order necessary to ensure that no contributions of order n are produced by the 
difference. Details of this calculation are described in Appendix A. The result is 

h(n) = n [ im) + 0(n°) . (29) 

^mi, (l + q4fl(l,q) 

Here Kfn is a standard geometric factor defined by 

Ku = (4.)-/ 2 ^ , (30) 

while V\ means the polynomial 

P 1 ( (Z 4 ) = 4-m(l + (Z 4 ). (31) 
The integral 7(1, q), defined by 



(to) r (m) 



p J q 



1 



IM= , ,2 , 747-, o . ,,, . -7T, ■ (32) 



{p' Z +r){\p' + p\ 2 + W + q\ 4 ) 



is the counterpart of Ma's (1973) critical "elementary bubble" n(fc) for the present 
case of m-axial LPs. 

The function I(p,q) satisfies scaling relations similar to those of G^(p, q) in 
equation (14): 

I(p, q) = p- £ 7(1, qp- 1 ' 2 ) = q- 2i I (pq' 2 , l) . (33) 

Further properties of this function are discussed in Appendix B. 

The O(n) term of 7 2 (n) can be worked out in a similar fashion. To calculate it, 
we must add and subtract the Taylor expansion of the corresponding inner integral 
to fourth order in one of its variables (see Appendix A). Proceeding in this manner 
yields 

Krn V2 (P 2 ) 

2 (rfM + 40( 1 ))to(to + 2) J P (1 +p 2 f J(p, 1) 
with the polynomial 

V 2 (p 2 ) = 3(8 - m)(6 - to) + 5(m 2 + 2to - 96)p 2 (35) 
+ (to 2 + 50to + 144)/ - m(m + 2)p 6 . 



j a(n) = 2 __ ^ / !™J + (n°) (34) 
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Using the expressions (29) and (34) for the integrals I\{n) and l2(n), we now 
solve the equations (27) and (28) for r^l and This gives 

1 2K rll /•(") V^) 1 _ 2 

»£2 = — / 7— — ttt 77: — r + 0(n ), 36 

n m J q (1 + g 4 ) 3 1(1, q) 

1 VL2 I K m /■(*) P 2 (p 2 ) 



.njW l_ + (n- 2 ) (37) 

1 (l+^/( P ,l) + ° (n j - 



2 4 n Am(m + 2) 

The implied 1/n expansion of 77^4 follows with the aid of relation (19) for the 
1/n coefficient rfl\ ; it reads 



= 1 2K m [W V 2 (p 2 ) 1 , _~ 

?7L4 nm(m + 2)J p (l+p^f I(p,l) +U[U h 



(38) 



In equations (37) and (38) we have expressed the 1/n coefficients in terms 
of an integral over the perpendicular momentum p. Upon exploiting the scaling 
property (33) of the integral I(p,q), we could rewrite this as an integral over the 
parallel momentum q. The transformed integrand would contain the polynomial 
?M? 4 ) = q 12 T y 2(q~ 4L )- Likewise, the (/-integral in equation (36) can be recast as 
an integral over p. 

The above equations (36), (37) and (38), giving the expansions of the correlation 
exponents 77L2, 7]L4 and the anisotropy exponent 9 to order 1/n for general values of 
(m, d) with 2 + 77i/2 < d < 4 + m/2, are the main results of this paper. In the following 
we further analyse these results for special choices of (m,d). Figure 1 illustrates the 
region between the upper and lower critical lines in which nonclassical behaviour is 
expected and displays the special cases to be discussed below as full circles. 



4. The limiting cases of critical points and isotropic Lifshitz points 

In this section we check that our results for general m, equations (36)-(38), comply 
with known series expansion to 0(l/n) for the cases m = of CPs (Ma 1973, Abe 
and Hikami 1973) and m = d of isotropic LPs (Hornrcich et al 1975b). In either 
case we must take limits D — > of Z?-dimcnsional integrals to determine the 0(l/n) 
coefficients, where D = m or fh. To do this, we use the relation 



Yim^Jd D xg(x;D) = hrn S D jdx x°- x g(x; D) = 3 (0;0) . 



(39) 



Consider first the CP case m —> 0. Here only the correlation exponent t]l2 is 
physically meaningful. From equation (36) one easily derives 

-l^im 1 ^ ^' (40 » 

where Jd(l, 1) is a standard one-loop integral, corresponding to a special case of the 
quantity Jo(f, t) defined by equation (A.l) of Appendix A. Its explicit value is given 
in equation (A. 9). Inserting it into the last equation we immediately recover the 
familiar result for 77 first derived by Ma (1973) and Abe and Hikami (1973): 

4(4-d)T(d-2) 1 _/„ , 

71 55 ""1"=° = dT{2 d/2) THd/2 - 1) fp n + ^ ) ' (41) 
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Figure 1. Lines of upper and lower critical dimensions d*(m) = 4 + m/2 and 
d((m) = 2+m/2, respectively. In the shaded region, bounded by portions of these 
lines together with the condition < m < d, nonclassical critical behaviour is 
expected to occur. The full circles indicate the special cases discussed in section 7. 
Also shown are the lines d = 3 + m and d = 1 + m on which the free propagator 
simplifies to a simple exponential and incomplete gamma function, respectively 
(see section 7). 



The limit m — > d, corresponding to an isotropic LP, can be handled in a similar 
fashion. Now only the correlation exponent T)L4 is physically significant. Equation (38) 
yields 

2K d .. 7> 2 (0) 1 



lim n Li 

m—>d 



d(d+2) Jr d mi)n +0 ^ 



(d-8)(d-6) 



Kd - + 0(l/n 2 ) 



(42) 



d(d + 2) J d (2,2)n 

with Jd(2, 2) given by equation (A. 10). Upon making several transformations of 
Gamma functions, we recover the result of Hornreich et al (1975b) for 77^4, namely 

. 2 sin(7rd/2)r[(d-3)/2] 1 



3(8 - d) 2° 



7T 3 / 2 d(d + 2) r(d/2) n 



+ 0(l/n 2 ). (43) 



5. Consistency with e expansion about the upper critical dimension 

In (Dichl and Shpot 2000, Shpot and Diehl 2001) the e expansions about the upper 
critical dimension d* = 4 + m/2 of all critical, crossover, and the usual correction- 
to-scaling exponents have been obtained for general m and n. These series can be 
expanded in powers of 1/n to 0(l/n). Conversely, considering the limit of small e, 
our above 0(l/n) results for general (m,d) can be expanded in powers of e. The 
resulting two double series expansions in e and 1/n of each exponent should agree. 

We shall work out the 0(e 2 /n) contributions implied by the 1/n expansions 
(36) and (38) for general m, expressing them in terms of multi-dimensional integrals. 
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Unfortunately, these integrals are in general rather complicated and not necessarily 
analytically tractable. For this reason, we will content ourselves here with verifying 
the consistency of the series expansions of rj L2 and rj Li in e and 1/n for the special 
choice m — 2. Owing to the simple form which the scaling functions of the position- 
space propagators take in this case both for n = oo and the Gaussian theory (simple 
exponentials, see Mergulhao and Carnciro (1999), Diehl and Shpot (2000), Shpot 
and Diehl (2001)), the 0(e 2 /n) terms of the series expansions can be worked out 
analytically in a straightforward fashion (see below). 

One source of e dependence in equations (36) and (38) is the function I(p, q). This 
is just the one-loop Feynman integral associated with the four-point graph XIZX 01 
the usual <p 4 theory for m-axial LPs. Its Laurent expansion in e reads 

I(p, q) = ^ c_r [1 - 6 F(q 2 /p)} + O(e) , (44) 

where c_i and minus the scaling function F represent the residuum and finite part of 
1(1, q), respectively. 

In the case of r]L2, a second source of e dependence is the function Pi(q 4 ) (sec 
equation (31)). Upon decomposing it as 

Vi(q 4 )=V[ 0) (q 4 ) + er{ 1) (q i ) (45) 

with 

7fV) = y(l + 9 4 )-4<z 4 (46) 

and 

V i 1 1 \q 4 ) = l + q 4 , (47) 
we arrive at the e expansion 

„ L2 . i ^ r ^V)±iffv> tl + « F(rt] + (f!, i) . (48) 

n m c-i J q (1 + q 4 ) 6 1 J \ n n z J 

The result suggests the presence of a term linear in e. However, the e expansions 
of the correlation exponents 77^2 and 77^4 are known to start at order e 2 , so this term 
must vanish. It does indeed, since 

r°° T>^(n 4 ) 

<TT7P = °- <49 » 

Hence we have 

n e> 4iT 4 _ m/2 f^ r[ a \q 4 )F(q 2 ) + r[ 1 \q 4 ) /e 3 £ 2 X 

VL2 -n (8-m)c.J q (TT^j3 + U'^J- (50) 

Turning to the expansion (38) of 77^4, we note that the polynomial "P 2 (p 2 ) in its 
integrand docs not depend on e. However, the integration measure contains a factor 
p~ e which cancels the one originating from I(p, 1) (see equation (44)). It follows that 

1 2K m e /■<*') T 2 ( P 2 ) 



VL4 = - . " — / n? 2 5 [l + eF(l/p)] + -,- ■ (51) 
nm(m + 2)c-i J p {l+p z )° \n n / 

Similarly as above, the contribution from the 1 in the square brackets vanishes because 

W- m/ 2^!) =0 (52) 

(1 +/J 2 ) 



/' 

Jo 
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for general m. Thus the contribution linear in e is zero, and the e expansion of the 
large-n result (38) becomes 

2if m r^ V 2 (pi)F(l/p) /e 3 e 2 x 

7 ? L4 = ? ro\ / — TT, — 2V5 hOI — , — )■ 53 

n m{m + 2)c_ij p (l+/r)° Vn n/V 

Equations (50) and (53) give the contributions of 0(e 2 /n) implied by our large-n 
results for general m. They should agree with those obtained from the e-expansion 
results of (Diehl and Shpot 2000, Shpot and Dichl 2001). Let us verify this explicitly 
for the biaxial case m = 2, for which both the residue c_i as well as the function F 
appearing in the Laurent expansion (44) may be gleaned from Mergulhao and Carneiro 
(1999). Taking into account that the expansion parameter e|| , called e by these authors, 
corresponds to = 2e in our notation,* one sees that 

c-iU 2 =3^2, (54) 

F (« 2 )L=2 = y arctan ( 2 /'? 2 ) + \ ln C 1 + ? 4 /4) - 1 - In 4 . (55) 

Upon substituting this into equations (50) and (53), one can perform the required 
integrals to obtain 

n i2( ™ = 2) = ^ + 0(^), (56) 

^( m = 2)=-^ + 0(£,£). (57) 

Expanding the correct 0(e 2 ) results (Sak and Grest 1978, Mergulhao and Carneiro 
1999, Shpot and Diehl 2001) to order 1/n yields identical results. 

Since the Laurent expansion of the integral I(p, q) to order e° is also explicitly 
known when m — 6 (where d* = 7) (Mergulhao and Carneiro 1999), a similar 
consistency check should also be possible in this case via analytical, albeit somewhat 
more complicated, calculations. 

A complete proof of agreement of the results of the 1/n and e expansions for 
general m lies beyond the scope of the present work. The main reason why we have so 
far been unable to generalize the foregoing consistency check to general values of m is 
our lack of knowledge of a sufficiently simple closed form for the finite part of I(p, q). 



6. Consistency with expansions about the lower critical dimension 

The 1/n expansions (36) and (38) should hold down to the lower critical dimensionality 
de, which is given by the line de(m) = 2 + m/2 in our case (unless the LP gets 
destroyed by fluctuations). Considering dimensions d = df + eg slightly above de, we 
can additionally expand in ee- Conversely, any pertinent series expansion in eg that 
is available can be expanded in powers of 1/n. The resulting pairs of double power 
series in 1/n and e must agree. 

Series expansion about the lower critical dimension can be obtained for systems of 
the kind we are concerned with here — namely, systems having a low-temperature phase 
with spontaneously broken 0(n) symmetry — by analysing appropriate nonlinear sigma 

* Note that there is a misprint in the corresponding formula of Shpot and Dichl (2001) which precedes 
its equation (69). The correct relation between and e = d*(m) — d is as given above in the main 
text. 
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models (Polyakov 1975, Brezin and Zinn-Justin 1976a, 1976b, Bardecn et al 1976). 
There are two sorts of results we can compare with: q expansions for the CP case 
m = (Brezin and Zinn-Justin 1976b) on the one hand, and for m-axial LPs (Grest 
and Sak 1978), on the other hand. By investigating an appropriate generalization of 
the conventional nonlinear sigma model, the latter authors produced one-loop results 
for the three special values to = 1, to = 2, and m = 4. 

To show that our 0(l/n) results (36)-(38) are compatible with these ee 
expansions, we proceed in much the same way as in section 5. The integrals I(j>, 1) 
and 1(1, q) in the denominators of these equations have a pole at ee = 0. According 
to equations (B.8) and (B.23), 

I(p, 1) = 1(1, q) + 0(e?) = (4 7 r)- 1 -W4 £|W| 1 + o(e °) . (58) 

Substituting this into equations (36)-(38), one can perform the remaining integrals. 
Unlike the case of small e, nonzero contributions linear in ee appear. We obtain 

VL2 = ^- + 0[^,— =- + 0(— ,— (59) 

and 

These results are consistent both with Grest and Sak (1978) as well as with Brezin 
and Zinn-Justin (1976b). Since they were derived for general to, the above critical 
exponents are independent of to to the indicated order in 1/n and e. Note that the 
0(et/n) correction to the anisotropy exponent vanishes. 

From the structure of the ee expansions it is clear that the O(q) terms of r]L2 and 
?7i4 must be inversely proportional to n — 2. Hence we can generalize Grest and Sak's 
(1978) results for m = 1, 2 and 4 to conclude that 

VL2(m) - \vL4(m) + 0(e\) = + O(ej) , (61) 

e(m)= 1 - + 0(ej), (62) 

for general m. 

As one sees, the series coefficients of the terms linear in eg arc independent of m. 
This is a feature they have in common with the coefficients of the 0(e) contributions 
of all bulk critical exponents of the m-axial LPs we are concerned with here (Shpot 
and Diehl 2001). For the e expansion of the above exponents this is trivially fulfilled 
since their 0(e) terms are zero. However, other exponents such as those associated 
with the parallel and perpendicular correlation lengths do have O(e) contributions, 
whose coefficients do not depend on m. 



7. Special cases 

In this section we further exploit our results for general m, equations (36), (37) and 
(38), by considering special cases. Of primary physical interest clearly is the case 
(to, d) = (1, 3) of a three-dimensional system with uniaxial anisotropy. Unfortunately, 
this case appears too difficult to allow a completely analytic calculation of the 0(l/n) 
terms. Let us therefore first consider some simpler cases, before returning to it. 
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7.1. The case m = 1, d = 4 

For this choice of (to, d), the required calculations can be done analytically to obtain 
closed expression for the 0(l/n) coefficients of the exponents r)L2, t\l\ and 6. 

To see this, recall that the free propagator in position space takes a simple 
exponential form on the whole line d = 3 + m (Mergulhao and Carneiro 1999, Shpot 
and Dichl 2001). On it, the scaling function & m ,d(v) of the Fourier back transform 

G^(r,z)=r- 2+ ^ m , d (v), v^z/V^, (63) 

of the free momentum-space propagator G^\p, q) introduced in equation (16) 
simplifies to 

$ m ,3+m(v) = (4^)- 2+e e-^ 4 , (64) 

whilst away from it, it generally is a difference of two generalized hypergeometric 
functions, a so-called Fox-Wright i<J/i function (Shpot and Diehl 2001). 

Previously, this simplifying feature was exploited in the context of the e expansion 
in two different ways: Mergulhao and Carneiro (1999) fixed the codimension m = d—m 
at TO = 3 and took m = 2 — 2e (and hence d = 5 — 2e) to expand about the point 
(to, d) = (2, 5). Shpot and Diehl (2001) performed a two-loop RG analysis for general 
fixed to and dimensions d = d*(m) — e. Owing to the simple form (64), both types 
of calculations could be performed analytically for to = 2, as well as for m = 6 where 
one can benefit from similar simplifications at d* (6) = 7. Despite the difference of 
the two procedures, they yielded consistent results for the e expansions of the critical 
exponents about (to, d) = (2,5) to 0(e 2 ).jj 

Since the e expansion is asymptotic, it gives us direct information only about the 
behaviour in the immediate vicinity of the line of upper critical dimensions d*(m). In 
order to derive from it reliable predictions for the values of critical exponents in d = 3 
dimensions, one must extrapolate, for example, in the case m = 1 of a uniaxial LP 
down to e = 3/2. This is not normally possible unless sophisticated extrapolation and 
resummation methods are employed. The 1/n expansion, on the other hand, does not 
require e to be small and hence enables us to move further away from d* . Here we 
exploit the simple form (64) to gain information about the behaviour at (to, d) = (1, 4), 
corresponding to a distance of e = 1/2 from d*(l) = 9/2. 

It is no complicated matter to compute the function 7(1, q) for to = 1 and d = 4; 
one obtains 



which yields I(p, 1) = p^ 1 ^ 2 7(1, p -1 / 2 ) by virtue of the scaling property (33). Upon 
inserting these expressions into equations (36)-(38), one can perform the required 
integrals in a straightforward manner. This gives the following results for the 
1/n coefficients of the exponents n L2 , and n Li : 

^ = ^-0-1021, 

0(1) =-ir77^-°- O272 < \ = 4. ((»(») 

ft Verifying the consistency of Mergulhao and Carneiro's (1999) series expansions to 0(e 2 ) about 
(m, d) = (2, 5) and (6, 7) with those of Shpot and Diehl (2001) requires the correction of two minor 
mistakes in the former reference; see section 4.2 of Shpot and Dichl (2001) for details. 
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Despite the smallness of the 0(l/n) corrections, these results provide clear 
evidence for the fact that the critical exponents 77^,2, VL4 an d 9 have nonclassical 
values below the upper critical dimension. This conclusion is in full accord with 
previous work based on the e expansion (Dichl and Shpot 2000, Shpot and Diehl 
2001). 

7.2. The special case m = 4, d = 5 

There is another line on which the scaling function <3? mj d simplifies: For d = m + 1 
it reduces to an incomplete gamma function (see equation (18) of Shpot and Diehl 
(2001)). Specifically for d = m + 1 = 5, it becomes the elementary function 



4> 



Its simplicity enables us to compute the integral 1(1, q) for this choice of m and 
d without much difficulty. We obtain 



id=5 



1 



1 



In 



1 + g 4 /4 



(68) 



.9 2 (9 4 + 3). 

and a corresponding result for I(p, 1) = p _1 /(l,p -1 / 2 ). These expressions can be 
substituted into equations (36)-(38) and the required integrals then evaluated by 
numerical means. This yields the values 



0(1) 



0.314 , 
-0.0728. 
0.045 , 



> m = 4 , d ■■ 



(69) 



for the 1/n coefficients of r)L2, 9 and t]la- Remarkably, the 1/n coefficient of 77^4 is no 
longer negative, as it was both for infinitesimally small e and at (m,d) = (1,4). 



7.3. The case m = 1, d = 3 

We now turn to the uniaxial, three-dimensional case (m,d) — (1,3). Unfortunately, 
we have not been able to evaluate fully analytically the required integrals of the 
1/n coefficients. 

Let us start from equation (32). The (rh = 2)-dimensional integration over the 
perpendicular momentum p' can be performed in a straightforward fashion, giving 

- 2q'\q' + qf + A + [g' 4 + (q< + qf + l] A 1 / 2 " 



i<i -* ) =^L dq ' A ~ 1 "' n 



where A stands for the expression 

A(q', q)=[l + (q' + qf] ' + (1 + 9 ' 4 ) 2 - 2q'\q' + qf - 1 . (71) 

The integration over q' can be regarded as an integral in the complex q' plane 
and after some work be shown to reduce toff 

dq' 



1° ^/(l + 4g' 2 g 2 ) [1 + \(4q> 2 + q 
ft We are indebted to S. Rutkcvich for this calculation. 



(72) 
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This in turn can be expressed as a complete elliptic integral or Gauss hypergeometric 
function to obtain 

-V4/1 , „4x-l/2 „ n l.i.i.2\ 



1(4 + ^-^(1 + , 4 )- 1/2 



2 ^i(i^;i;fc 2 ) 



with 



e = 



i 



9 



(73) 



(74) 



(l + ^)(l + 4 /g4) 1 /2_ 

Known quadratic and linear transformation formulae for the hypergeometric 
function 2-F1, such as equations (15.30.30) and (15.3.4) of Abramowitz and Stegun 
(1972), enable us to cast the above integral in the following two equivalent forms 

/(M) = ^^^(1,1;^), {l + qi * {A + qiy (75) 

= ^ SFl( i,3 ;1; _ w); . = ^_. (76) 

Upon substitution of either one of them into our general expressions (36)-(38) for the 
0(l/n) coefficients (along with their counterparts for I(p, 1) = p _3 / 2 /(l,p -1 / 2 )), the 
remaining integral over q can easily be performed numerically. The results are 



VL2(m = 1, d = 3) ~ 0.306 - + 0(n~ 2 ) , 

n 



6{m = l,d= 3) 



1 -0.0487- + 0(n~ 2 ) 

2 n 



1 



r?L4(m = 1, d = 3) ~ 0.223 - + 0(n~ 2 ) . 



(77) 
(78) 
(79) 



Just as in the previous case of m = 4, d = 5, and unlike the coefficient of the 
0(e 2 /n) term, rf£l is positive. This suggests a tendency of 77^4 to change from small 
negative values of d < d* to positive ones as d is lowered further. 



8. A universal scaling function and related amplitude 

So far we have focused our attention on the calculation of critical exponents. However, 
the 1/n expansion can also be employed to compute other universal quantities such 
as universal amplitude ratios and scaling functions. As an example, we here compute 
the leading nontrivial term of a universal scaling function associated with the energy- 
density cumulant at the LP and a related amplitude. 

According to general scaling arguments and RG work (Diehl and Shpot 2000, 
Shpot and Diehl 2001, Diehl et al 2003a, 2003b), the leading singular part of this 
function is expected to take on large length-scales the asymptotic scaling form 

G^ip, q) = <</>V>p7L g « ^ P~* L/UL2 ^(E2 qp- 8 ) . (80) 
Here the quantity on the left-hand side is the Fourier transform of the cumulant 
(<j) 2 (x)(j) 2 (0)) cum , where "sing" means singular part. Further, aj, and 1/L2 are the 
usual critical indices of the specific-heat and perpendicular correlation length; they 
are related to the exponents used above via 

ollIvl2 = 2/^L2 -fh-m0 = 2A^, . (81) 

Furthermore, E\ and E2 denote nonuniversal metric factors. The former can be varied 
by changing the normalization of the order-parameter field </>, the latter by modifying 
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the scale in which parallel momenta are measured, i.e. by multiplying the term (<9 2 </>) 2 
of the Hamiltonian by a factor a ^ 1. 

The scaling function \t mj£ ; should be universal, but depends of course on n. We 
choose its normalization such that ^> m ,d{0) = 1- To fix the scale of its argument, we 
could require that the logarithmic derivative of ^ m ,d{<^) at a reference value qo (for 
instance, qo = 0) takes a certain value. We do this simply by making the choice £^2 = 1- 
The 1/n expansion of ^f m ,d then starts at order (l/n)°. Taking into account that the 
function G^(p,q) is trivially related to G$(jp,q) (see e.g. chapter 29 of Zinn- Justin 
(1996)) and recalling equation (21), one sees that 

* TOld (q) = ^jj+0(l/n). (82) 

If we take the limit p — > at fixed q, the scaling form (80) must yield a p-independent 
contribution ~ q- a ^l v ^ . Hence the scaling function must have the asymptotic 
behaviour 

tf m , d (q - oo) « C(m, d) q-<W(""*) . (83) 

From equation (82) we conclude that the universal C(m, d) coefficient, to leading 
order in 1/n, is given by 

C{m,d) =C (0) (m,rf) + O(l/n) with C (0) (m, d) = . (84) 

7(1, 0) 

In Appendix B we show that the required two integrals can both be calculated 
explicitly for general values of m and d. The results can be found in equations (B.8) 
and (B.23). We refrain from giving the resulting rather lengthy expression for C(m, d) 
here, restricting ourselves to a discussion of some special cases of interest. 
Consider, first, the uniaxial case m = 1. Here our result becomes 

_ 2- d + 3 / 2 y^r [(2d - 3)/4] sin [(2 d - 5) tt/4)] [2 d + 16 sin(d7r/2)] 

u ' j_ cos(d7r)r(i/4)r[(d-i)/2] • [ ' 

As shown in figure 2, this coefficient is a smooth and finite function of d in the whole 
range de = 5/2 < d < d* = 9/2. It has the e expansion 

C (0) (l,9/2-e) = l + (7r-4-21n2)^ + 0(e 2 ). (86) 

The corresponding results for the biaxial case m — 2 are 

C (0) (2,rf) ^- 3/2 cot(^/2) 



r[(5 - d)/2] T[(d - 4)/2] - 4 5 -V 4) 



(87) 



T(d-7/2)_ 
and 

C (0) (2, 5 - e) = 1 + (In 2 - 1) e + 0(e 2 ) . (88) 

Again, this coefficient behaves smoothly between the upper and lower critical 
dimensions dg = 3 and d* — 5 (sec figure 2). 
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0.95 



Figure 2. Universal amplitude C(m, d) as a function of d for m = 1 (full line) 
and m = 2 (broken line) 



9. Summary and Discussion 

In this paper we have shown how to utilize the 1 /n expansion for the study of critical 
behaviour at m-axial LPs. We have been able to determine the 0(l/n) corrections of 
the correlation exponents 77^2 and 77^4, and of the related anisotropy exponent 9, for 
general values of to and d with < to < d and 2 + m/2 < d < 4 + to/2. What makes 
such calculations a challenge is a combination of two problems: the anisotropic scale 
invariance one encounters already at the level of the free theory, and the complicated 
forms of the propagators' scaling functions at the LP. 

To cope with these difficulties, it proved useful to employ a properly adjusted 
and generalized modification of the technique Vasiliev et al (1981a, 1981b) introduced 
to compute the critical exponents r] and v of the standard n-component |(/>| 4 model 
up to order n -2 . We believe that the present work may serve as a starting point 
for more ambitious studies based on the 1/n expansion. One question which might 
be systematically investigated in this manner is whether the predictions for scaling 
functions of anisotropic scale invariant systems made recently by Henkel (1997, 2002) 
have any significance in cases where fluctuation effects must not be ignored. 

Apart from yielding insights into the feasibility of the approach to such anisotropic 
scale invariant system, our work permits us to draw two interesting conclusions. First 
of all, it provides unequivocal evidence for the fact that the critical exponents r]L2, 
?7i4 and 9 differ for d < d* from their classical values. Although this is in complete 
accord with what the available dimensionality expansions about the upper and lower 
critical dimensions tell us (see sections 5 and 6), it goes considerably beyond these 
results because no extrapolation in d is involved in the 1/n expansion. This means 
in particular that the anisotropy exponent 9 of the three-dimensional uniaxial LPs 
should be nonclassical. 

Another remarkable aspect of our results is the interesting variation of the 0(l/n) 
coefficient of the exponent 77^4 with d. As illustrated in figure 3 for the uniaxial case 
to = 1, it first decreases to small negative values as d drops below d* , then becomes 
positive as d is lowered further, before it drops back to zero at the lower critical 



Large-n expansion. 



19 



€i expan 
sion 




0.3 O 



0.2 



0.1 



2.5 



3.0 



3.5 



O. 



"W 7 
4.0 



e expansion 



4.5 

e expansion 



Figure 3. Behaviour of the 1/n coefficients r7^ 2 (f , d) (open circles) and ?7^ 4 (1, d) 
(crosses) as functions of d. The 4 points displayed for either one of them 
correspond to results described in section 7. The thick lines near the upper and 
lower critical dimensions represent the limiting forms ~ e 2 and ~ eg, which the 
dimensionality expansions mentioned in sections 5 and 6 yield. The widely or 
densely dotted lines serve to guide the eyes. The exponent 77^(1, d) has a small 
negative value beneath the upper critical dimension d*(l) = 9/2 and appears to 
change sign somewhat below d = 4. 



dimension dg. If we accept the plausible hypothesis that the correlation exponents are 
continuous functions of d, then it seems reasonable to expect a qualitatively similar 
behaviour of the exponent 77^4. In other words, 77^4 should change sign somewhere 
below d*. In fact, there are other examples where the lowest nontrivial term of the 
dimensionality expansion about d* is negative, whereas extrapolations of higher-order 
calculations yield a positive value in d = 3 dimensions. This happens for instance 
in the case of the Ising model with quenched random bond disorder. Here the 
dimensionality expansion of the correlation exponent rj about d — 4 begins with a 
negative contribution ~ (\/4 — d) 2 (see e.g. Jayaprakash and Katz 1977), whereas the 
extrapolation of higher-order calculations yields a positive value for this exponent at 
d = 3 (Pclissctto and Vicari 2000). 

Let us close with a few cautionary remarks. Our results to 0(l/n) are 
mathematically well defined in the whole range where 2 + m/2 < d < 4 + to/2 and 
< to < d (the shaded region in figure 1). However, applying them to real three- 
dimensional systems requires some care. They should be combined with information 
from other sources such as dimensionality expansions, simulations, exact results and 
experiments, keeping in mind their limitations. For example, depending on whether 
to = 2 or to > 2, the lower critical dimension dp is dg = 3 or dt > 3, respectively, 
whenever n > 2. Hence, unless the 0(n) symmetry is explicitly broken, no m-axial 
LPs are expected to occur at nonzero temperatures in d = 3, which leaves us with 
the case of uniaxial LPs in three dimensions. On the other hand, contributions to the 
Hamiltonian that break the 0(n) symmetry are generically expected, for example, for 
magnetic crystals. An analysis of such spin anisotropies can be found in (Hornreich 
1979). 

For m > 2, one must also worry about another type of anisotropies: space 
anisotropies breaking the isotropy in the m-dimensional subspace. These give rise 
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to corresponding anisotropic terms quadratic in (d 2 (j)/dz a dz a i), which we also have 
not taken into account here. Their effects have been investigated recently within the 
framework of the e expansion (Diehl et al 2003b). 

Last, but not least, let us mention that the technique of matching asymptotic 
scaling forms, on which our above analysis was based, can be extended to determine 
the 0(l/n) corrections of the thermal exponents vl2, ^L4, 1l, and for m-axial LPs 
(Shpot et al 2005). 
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Appendix A. Large-n behaviour of momentum integrals 

In subsequent calculations we shall need the properties of the following D-dimensional 
momentum integral: 

r(D) 

J D (f,t)= k- 2 f\k + l\- 2t = V(f,t;D-f-t), (A.l) 

Jk 

where the function V is defined by 

V(c 1 ,c 2 ;c 3 ) = (A7r)- D ^a(c 1 )a(c 2 )a(c 3 ), a(a) = T ^ 7 ° l) , i = 1,2,3. (A.2) 

Let us consider the parameters / and t of the integral (A.l) as regular functions 
of some small parameter k. Choosing I > to be integer, we set 

f + t=j-l + K . (A.3) 

Then the integral Jd(/, t) is singular in the region of large k for small values of k. Its 
singular part is given by 

where K D was defined in equation (30), and (• • ■); are Pochhammer symbols: 

(w)i = T{ Z\ l) =w(w + l)---(w + l-l), Ho = l. (A.5) 
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Specifically for I = 0, I = 1, I = 2 and n — > 0, we have 

V sing (/,i ;J D/2- K ) = ^, (A.6) 



V sing (/, t; D/2 + 1 - k) = / *U= , (A.7) 



V sing (/, i; D/2 + 2 - «) = D( £l 2)R /(/ + 1) *(* + l)U=o • (A.8) 
Two relevant examples of the integral Jo(f,t) are 

1} " /J* ?sW = (4 ' r " 2 '-W^W T2{d ' 2 - 1] (A ' 9) 

and 

2) = t" rTT^rru - (*r)-* r> (< - 2 V ,A,0, 



« 4 (<r + l) 4 r(d-4) 

Now let us consider the momentum integrals I\ (n) and I 2 (n) of equations (27) and 
(28). We have to isolate their pole contributions at n = oo. For further convenience 
we write them as 

/■(m) iA^-m /-(m) 20 . i 

h(n)= -Z ^/ Gi( , - „. g 2 ) — r, (A.ll) 

n > J P \p + l\^ J q 'Vlp+ll 20 q J F(l,q)' V j 

/■(to) AA^/6-m ,(m) 2/0 j 

where Gi(Q 2 ) and G 2 (P 2 ) denote the scaling functions G^ as) (l,Q) and G^ as) (P, 1). 

Adding and subtracting the asymptotic forms of the arguments of Gi and G2 for 
p — > 00 and q — > 00 we write these functions as 

G 1 (q 2 p 2e \p + ir 29 ) ee Gl (q 2 + , 2 a p ) = X) i q 2s a p + R t (p, q), (A.13) 

s=0 S - a[q > 

1 d s G 2 (p 2 



G 2 (p 2 g 2 / e | Q + ir 2 /^G 2 ( P 2 + P 2 /3 q )=^4^^ip 2 ^ + i? 2 fe (Z ). (A.14) 
The deviations 

ap ^p 2e |p + ir 20 -l and P q ^q 2 l e \q + l\- 2 / e -1 (A.15) 

are of order p~ x and g _1 for large p and g. Including their successive powers enhances 
the ultraviolet convergence of p and q integrations in equations (A.ll) and (A. 12). 
The remainders R\ (p, q) and R 2 (p, q) of the Taylor expansions contain sufficiently 
high powers of a p and (3 q to make the integrations finite as n — > 00. These convergent 
contributions would be required for a calculation to 0(n~ 2 ). All we need now is 

W = E^*3,/ ^^^)+0(A (A.16) 

I a (») = f 4 Hm ^1 -1- + O(n0) . (A.17) 

y ' f^s\ ™j p * d( P 2 )* F(p,l) v ; v ; 
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The coefficients Jf and Jf are given by the integrals 

Js = / p 4A *- rn \p + lp 2A * a s p , s = 0,l,2, (A.18) 



r-(m) 

Jf=/ + * = 0,...,4, (A.19) 



which have simple poles in 1/n at n — oo. This fact allows us to take the limit n — > oo 
in the remaining integrals over q and p in equations (A. 16) and (A. 17). 

Making binomial expansions of the integer powers a p and /?* in equations (A.18) 
and (A.19), we can express Jf and jf in terms of the integral (A.l): 

Js = E(-l) s ^^ V f | - 2A - j0, A, + j0; | + 1 - , (A.20) 



2 " ' ' v J ' 2 2n 
-(_-2A^j,— ^— ,- + 2 



(A.21) 



where C| are the corresponding binomial coefficients. 

Comparing the arguments of the function V in equations (A.20) and (A.21) with 
their counterparts in equation (A. 4), we read off Z = 1 and k = r/^/C^n) for Jf, and 
/ = 2 with k = {ri ( H + 46>( 1 ))/n for jf . With the aid of the results (A.7) and (A.8) we 
then obtain the required leading large-n behaviour of Jf and Jf: 

Jf = ^^A s + 0(n°), (A.22) 

Jf = Km Bs + 0(n°) . (A.23) 

m(m + 2) „W + 4^(1) 

Here A s (with s = 0, 1, 2) and (with s = 0, . . . , 4) denote the sums 

^ = £(-ir j+1 Ci ("-2-0(1 + 0, (A.24) 

s s = E(-i)^ c i (y - 4 - 2 i) (y - 3 - 2 i) ( 2 + W + 2 ^') • (A.25) 

3=0 

A simple calculation gives: 

A = ^{4-d + m), A 1 = ^(7-d + m), A z = \> 

B = ^(8-to)(6-to), Bi = ^(16 - m)(66 - 7m) , (A.26) 
B 2 = 2(612- 58m + m 2 ), B 3 = 48(24 -m), B 4 = 384 . 

In conjunction with the relations (A.22) and (A.23), these results give us the 
coefficients Jf and Jf in equations (A. 16) and (A. 17). 

The limit n — > oo in these equations reduces the scaling functions in the integrands 
to those of the free theory. That is, we may use the large-n result 

rid) 1 

lim F(p, q) = I(p, q) = / (A.27) 
™^°° Jfc' (p' + g' ) (\p' + p\ 2 + \q' + qr) 
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for both scaling functions F(l, q) and F(p, 1) introduced in equation (24). 

Similarly, the scaling functions G\{q 2 ) and G 2 (p 2 ) in equations (A. 16) and (A. 17) 
reduce to 

Gf\q 2 ) = T ^ ¥ and G<°V) = (A.28) 

in this limit. 

We can now evaluate the sums in equations (A. 16) and (A. 17). Apart from overall 
factors, we have 

^ 8\ q d(q 2 Y 2(l + g 4)3 ' y A -^> 

f^f&.M. (A30) 
j^s\ P d{p 2 Y 2(l+p 2 )5 ' [ ' 

where the polynomials T > i(q 4: ) and T l 2(p 2 ) are given explicitly in equations (31) and 
(35) of the main text. Finally, equations (A. 16) and (A. 17), along with (A.22)-(A.23) 
and (A.29)-(A.30), yield the leading-order expressions of h(n) and hin) given in 
equations (29) and (34). 

Appendix B. The one-loop integral 

The one-loop Feynman integral I(j>, q) is an important ingredient of our calculations. 
Unfortunately, we did not succeed to evaluate it for generic dimensions and external 
momenta. Various simplified expressions which result for special values of m and d are 
presented in the main text. Below, we consider the simplifying cases of zero external 
momenta where closed analytic expressions are obtained for arbitrary m and d. 

Appendix B.l. The integral l(p,0) 
Consider the integral 

Am) Am) ^ ^ 

I{p ' 0) = J P , l ^WT? (R1) 

It is convenient to convert this integral to the coordinate representation where it is 
given by a Fourier transform of the squared free propagator: 

l(p,0) = J d"V J d m z [Gj^r,*)] e~ ipr . (B.2) 

Using the scaling representation (63) for G^ (r, z) and changing the integration 
variable z via z = Vy/r we obtain 

I(p, 0) = J d"V r - 4 +™/ 2 + 2e e - ipr J d m v $ 2 (u) . (B.3) 

The r integration is standard, giving 

d fn r r -A +m /2+2e -ipr = ^ ^2 £(e/2j 

r(2-m/4-e) F v ; 
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To calculate the integral over v, we use the Fourier representation (Diehl and Shpot 
2000, Shpot and Diehl 2001) 

Am) 

*(v) = (2^)-W2 / q r-n-2 K ml ^_ x {<?) e^ . (B.5) 

Jq 

This leads us to 

d m v $ 2 ( W ) = (27r)-*y q^Kl^tf ) . (B.6) 
The last integral is known for arbitrary m and d. We obtain 

/ d-t; * 2 (,) = (2.)"-^ 1 2 -^(l - e/2) ^| • (B.7) 
Multiplying both contributions yields 

^» = <^1r W2 >M^ P -.. ( B. 8) 

Appendix B.2. The integral 1(0, q) 

Let us consider another special case of the one-loop integral I(p, q), 

'(0,1)=/ / ■ 7Z I Tvl • (B.9) 



p 2 + ? 4 p 2 + (q + l) 4 ' 

Its treatment is somewhat more involved. We shall use some tricks employed by 
Mergulhao and Carneiro (1999) in a similar calculation. Making a partial fraction 
expansion 

1 



(a + b)(a + c) c — b \a + b a + cj 
we represent (B.9) as 

Integrating over p gives 

f(m) | 



(B.10) 



(B.ll) 



p 2 +q 4 p 2 + ( q+1 )4_ 

1(0, 1) = (^r^T(-u) J (g + 1)4 _ g4 {q Av -\q + 1\ 4 ") (B.12) 
with 

to m e ,„ 

i/= 1 = 1 . (B.13 

2 4 2 y ' 

It is useful to represent the difference in the last brackets via the elementary integral 

l 

a a - b a = -a (b-a) [ t f X ^ (B.14) 

v 'J [a + x(b - a)] 1 -" v ; 

o 

with a = q 2 , b = (q + 1) 2 , and a = 2u. The factor (& — a) cancels with a corresponding- 
term of the denominator (see equation (B.12)). We obtain 

m i ) = tpW — — -i _ f dx I _ (B.15) 

1 7 (4tt)™/2 J g q 2 + ( q + l)2j (q2 + 2qlx + x) 1 - 2v 1 ' 
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Here, we represent the momentum-dependent factors through Laplace integrals 



This yields 
1(0,1) = 



r( 



{-2v)T{-v) 



oo 

(a) J 



(B.16) 



(4 7r )m/2 r(! _ 2z y) 

1 oo oo 



X 



y dx J dy J dzz- 2v J {m \-^ 2+2 ^ +1 ^e^+ 2 'i lx +^ . (B.17) 





Performing the Gaussian integration over g we get 

1 oo oo 



7(0,1) = 



(47r) d / 2 r(l - 2i/) 



dx dy dz 



-1v 



2y + z 



(B.18) 



(2y + z)" 1 / 2 

ooo 

The integrals over t and z can be decoupled by a rescaling y — zt oi the integration 
variable y. This gives 

1 oo oo 

(At:)-^T{1-2v)J dX J (2t+iy^ J dZZ 6 ' { ^ } 



The z integral is of the form (B.16). We obtain 

1 oo 



7(0,1) 



(-2^)r(-^)r(e) 

(47r) d / 2 r(l - 2v) J J (2t+l) 





f dx f - dt j- (t 2 + t + x- x 2 ) 

J J 2t+lW2-e V > 



(B.20) 



Inside of the last brackets we add and subtract 1/4. Next, we introduce the new 
integration variables y = 2t + 1 and s = 2x — 1. A little bit of algebra then yields 

{-2v)T{-v)T{e) 2 2 < 



7(0,1) 



(47r) d / 2 T(l - 2v) 4 



J ds j dy y- m ' 2+ ^y 2 - s 2 y\ (B.21) 



l 

Here the last integration produces a Gauss hypergeometric function: 

l 



7(0,1) = 



{-2v)T{-v)T{e) 2 



2e-l 



(4 7r )d/2 r(l - 2z/) m-2 + 2e 

-l 

Performing the remaining integral we finally obtain 

■r(e)r(l-e)2 



ds 2 F 1 (e,^^; 2 ^tM;s 2 ) . (B.22) 



7(0, q) = q- 2€ (4n)- d / 2 F( v) rwrn _ ^ 0-2+2, 



T(l - 2v) 



Hi 4 r[(2 + m + 2e)/4] 



_r(3/2-e) ' 2-m-2e r[(2 + m-2e)/4]J ' (R23) 

where we reintroduced the initially suppressed trivial dependence on q, taking into 
account that v = 1 — to/4 — e/2. 
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